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Abstract. Some Bianchi type IX viscous fluid cosmological models are investigated. 
To get a solution a supplementary condition between metric potentials is used. 
The viscosity coefficient of bulk viscous fluid is assumed to be a power function of 
mass density whereas the coefficient of shear viscosity is considered as proportional 
to scale of expansion in the model. The cosmological constant A is found to be 
positive and is a decreasing function of time which is supported by results from 
recent supernova observations. Some physical and geometric properties of the models 
are also discussed. 
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1. Introduction 

Bianchi type IX cosmological models are important and interesting in 
the sense that these models allow not only expansion but also rotation 
and shear and in general are anisotropic. Bianchi type IX universes 
include closed FRW models. The homogeneous and isotropic FRW 
cosmological models which are used to describe standard cosmological 
models, are particular case of Bianchi type I, V, and IX space-times 
according as the constant curvature of the physical three-space, t = 
constant, is zero, negative or positive. In these models, neutrino viscos- 
ity explains the large radiation entropy in the universe and the degree of 
isotropy of the cosmic background radiation. The standard cosmological 
models are too restrictive because of the insistence on the isotropy of 
the physical three-space and several attempts have been made to study 
non-standard cosmological models e. g. MacCallum (1979); Narlikar 
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and Kembhavi (1980); and Narlikar 1983). It is therefore interesting 
to carry out the detailed studies of gravitational fields which are de- 
scribed by space-time of various Bianchi types. Vaidya and Patel (1986) 
have studied spatially homogeneous space-time of Bianchi type IX and 
they have outlined general scheme for the derivation of exact solu- 
tions of Einstein's field equations in presence of perfect fluid and pure 
radiation fields. Krori te al. (1990); Chakraborty and Nandy (1992) 
have investigated cosmological models of Bianchi type II, VIII and 
IX. There are many other researchers (Uggla and Zur-Muhlem, 1990; 
Burd, Buric Ellis, 1990; King, 1991; Paternoga and Graham, 1996), who 
have studied Bianchi type IX space-time in different context. Recently 
Bali and Dave (2001, 2003) have investigated Bianchi type IX string 
cosmological models. 

Models with a dynamic cosmological term A(t) are becoming pop- 
ular as they solve the cosmological constant problem in a natural way. 
There is significant observational evidence for the detection of Ein- 
stein's cosmological constant, A or a component of material content of 
the universe that varies slowly with time and space and so acts like 
A. Recent cosmological observations by High -z Supernova Team and 
Supernova Cosmological Project (Garnavich et al, 1998; Perlmutter et 
al, 1997, 1998, 1999; Riess et al, 1998; Schmidt et al, 1998) suggest 
the existence of a positive cosmological constant A with magnitude 
A(G7t/c 3 ) ~ 10~ 123 . These observations on magnitude and red-shift 
of type la supernova suggest that our universe may be a accelerating 
with a large function of the cosmological density in the form of the 
cosmological A-term. Earlier researchers on this topic, are contained 
in Zeldovich (1968), Weinberg (1972), Dolgov (1983, 1990), Bertolami 
(1986), Ratra and Peebles (1988), Carroll, Press and Turner (1992). 
Some of the recent discussions on the cosmological constant "problem" 
and consequence on cosmology with a time-varying cosmological con- 
stant have been discussed by Dolgov (1993,1997), Tsagas and Maartens 
(2000), Sahni and Starobinsky (2000), Peebles (2002), Padmanabhan 
(2003), Vishwakarma (1999, 2000, 2001, 2002), and Pradhan et al 
(2001, 2002, 2003). This motivates us to study the cosmological models 
in which A varies with time. 

The majority of the studies in cosmology involve a perfect fluid. 
However, observed physical phenomena such as the large entropy per 
baryon and the remarkable degree of isotropy of the cosmic microwave 
background radiation suggests analysis of dissipative effects in cosmol- 
ogy. Furthermore, there are several processes which are expected to give 
rise to viscous effects. These are the decoupling of neutrinos during the 
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radiation era and the decoupling of radiation and matter during the 
recombination era. Bulk viscosity is associated with the GUT phase 
transition and string creation. The model studied by Murphy (1973) 
possessed an interesting feature in that the big bang type of singularity 
of infinite space-time curvature does not occur to be a finite past. 
However, the relationship assumed by Murphy between the viscosity 
coefficient and the matter density is not acceptable at large density. 
The effect of bulk viscosity on the cosmological evolution has been 
investigated by a number of authors in the framework of general theory 
of relativity (Padmanabhan and Chitre, 1987; Johri and Sudarshan, 
1988; Maartens, 1995; Zimdahl, 1996; Pradhan, Sarayakar and Bee- 
sham, 1997; Kalyani and Singh (1997; Singh, Beesham and Mbokazi, 
1998; Pradhan et al., 2001, 2002). This motivates to study cosmological 
bulk viscous fluid model. 

Recently Bali and Yadav (2002) has investigated Bianchi type IX vis- 
cous fluid cosmological models. Motivated by the situations discussed 
above, in this paper, we shall focus upon the problem with varying 
cosmological constant in presence of bulk and shear viscous fluid in 
an expanding universe. We do this by extending the work of Bali and 
Yadav (2002) by including varying cosmological constant and the coef- 
ficient of bulk viscosity as function of time. This paper is organized as 
follows. The metric and the field equations are presented in section 2. In 
section 3, we deal with the solution of the field equations in presence of 
viscous fluid. Section 4 includes the solution of some particular models 
whereas in section 5, we deal with some special models. In section 6, 
we have given the concluding remarks. 



2. The metric and field equations 

We consider the Bianchi type IX metric in the form 

ds 2 = -dt 2 + A 2 dx 2 + B 2 dy 2 + (B 2 sin 2 y + A 2 cos 2 y) dz 2 

—2A 2 cosy dx dz, (1) 
where A and B are functions of t only. 

The Einstein's field equations (in gravitational units c = 1, G = 1) 
read as 

R\ - \r£ + Agj = -8nT? (2) 

where R{ is the Ricci tensor; R = g %3 R, L j is the Ricci scalar; and T- is 
the stress energy-tensor in the presence of bulk stress given by Landau 
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and Lifshitz (1963) 

T i = (P + P)ViV j + pg\ - r] (v{. + v j -i + v j v l Vi-i + ViV l v j .^ 

2 



(3) 



Here p, p, r\ and £ are the energy density, isotropic pressure, coefficient 
of shear viscosity and bulk viscous coefficient respectively and v\ is the 
flow vector satisfying the relations 



9ijV v J 



(4) 



The semicolon (; ) indicates covariant differentiation. We choose the 
coordinates to be comoving, so that 



= v 2 



v 3 ,v 4 



(5) 



The Einstein's field equations (2) for the line element (1) has been set 
up as 



-8vr 



-8vr 



r 9 A * 


(<- 




V - 2ri— - 




\ 9 B ± 






p-2r]— - 







2B, 



B 

A 44 



B A 2 1 3A 2 , , 

+ ^ + ^"4^ + A ' (6) 



A4B4 U44 A 2 . . 



8vrp 



+ 



B. 2 



1 A 2 



(8) 



AS 5 2 B 2 AB 4 
where the suffix 4 at the symbols A and B denotes ordinary differ- 
entiation with respect to t and 6 is the scalar of expansion given 
by 

= «V (9) 



3. Solution of the field equations 

Equations (6) - (8) are three independent equations in seven unknowns 
A, B, p, p, rj, £ and A. For the complete determinacy of the system, we 
need four extra conditions. 

Firstly we assume a relation in metric potential as 

A = B m (10) 

and secondly we assume that the coefficient of shear viscosity is pro- 
portional to the scale of expansion, i. e., 

7?oc6» (11) 
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where m is a real number. 
Eqs. (6) and (7) lead to 

-B44 B4 2 A44 A 
~W + ~W ~ ~A AB~ ~B 4 'B 2 



A 2 1 ( A A B A 

—7 + -^7T = 167T77 — — 



V A B 



Condition (11) leads to 



T) = t 



A 4 2B 4 
~A + ~B 



where t is a proportionality constant. 

Equations (12) together with (10) and (13) leads to 



BB U + aB 4 2 = 



which can be rewritten as 



B 2(m-1) 



1 



(1 — m) (1 — m 



2B 



2m-3 



dB \J )~r B KJ ) { i_ m) (i- m )B 



where 



and 



a = (1 + m) - 16tt£ 



(to 2 + to — 2) 



(1 — to) 



Ba = f(B) 



From (15), we obtain 

dB\ 2 
~dt J 



B 2(m-1) 



+ 



P 



(1 - m)(m + a - 1) B 2a a(l - m 



(12) 



(13) 



(14) 

(15) 

(16) 
(17) 

(18) 



where (3 is a constant of integration. After a suitable transformation of 
coordinates, the metric (1) reduces to the form 



ds 2 = - 



ji2(m-l) 



+ 



P 



1 



(1-to) (m + a- 1) T 2a a(l - m) 



dT 2 +T 2m dx 2 + 



T 2 dy 2 + (T 2 sin 2 y + T 2m cos 2 y) dz 2 - 2T 2m cos y dx dy, (19) 
where B = T. 

The pressure and density for the model (19) are given by 



8ttb = K] T 2 (™- 2 ) + ^ + PKs 

onp ^ ^2 ^ ^2(a+l) 



+ 
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8vr£(m + 2) 



1 



j<2(m-2) 

(m + a- 1)(1 



+ 



m 



1 

y2(a+l) a (l _ m )T 2 



8vrp = K 4 T + T2(a+1) + ^ + A, 

where 

(-m 2 )(64rf+21)-m(64rf~3a+6)+(1287r^-3a+15) 
12(ra+a-l)(l-m) ' 
_ m 2 (167rl+3)+rn(167rl)-327rl 
2 3a(l— m) ' 

K 3 = i[(-m 2 )(16vr£ + 3) - m(167r^ - 3a) + (32vr^ + 3a)], 

m 2 +m(a+6) — (a— 5) 
4(m+a-l)(l-m) ' 



-A, 
(20) 
(21) 



K- 



K 4 



K _ (-m)(a+2)+(a-l) 
^5 ~ a(l-m) 

For the specification of £, we assume that the fluid obeys an equation 
of state of the form 



V = IP, 



(22) 



where 7(0 < 7 < 1) is constant. 

Thus, given £(t) we can solve for the cosmological parameters. In most 
of the investigation involving bulk viscosity is assumed to be a simple 
power function of the energy density (Pavon, 1991; Maartens, 1995; 
Zimdahl, 1996) 

m = (23) 

where £0 an d n are constants. If n = 1, Equation (23) may correspond 
to a radiative fluid (Weinberg, 1972). However, more realistic models 
(Santos, 1985) are based on n lying in the regime < n < \. 
On using (23) in (20), we obtain 



J<2 y2(a+l) 



8^0/0™ (m+2) 



y2(m-2) 



+ 



(m + a-l)(l-m) T 2 ( Q+1 ) a(l - m)T 2 



-A. 

(24) 



3.1. Model I: 



Solution for £ = £ 



When n = 0, Equation (23) reduces to £ = £0 = constant. Hence in 
this case Equation (24), with the use of (21) and (22), leads to 

»7r(l + 7)/? = (Ki + A4) T > + ^ + T 2( a+ i) 
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+87r£o(m + 2) 



1 



j<2(m-2) 



+ 







(m + a - 1)(1 - m) T 2 ( a+1 ) a(l - m)T 2 



(25) 

Eliminating p{t) between Equations (21) and (25), we have 
(1 + 7 )A = ( Kl - 7 K 4 ) T^-V + ^ ~J Kb) + HK "~ 7(2m + 1)} 



J>2 



+87re (m + 2) 



\ 



^2(ra-2) 



+ 



/3 



(m + a - 1)(1 - m) T 2 ( a+1 ) a(l - m)T 2 



(26) 



3.2. Model II: 



Solution for £ = £ p 



When n = 1, Equation (23) reduces to £ = £oP- Hence in this case 
Equation (24), with the use of (21) and (22), leads to 



8yrp 



l + 7-£o(™ + 2) 



y2(m-2) p 1 

(ra+a-l)(l-ra) ~r T 2 («+ 1 ) a(l-m)T 2 



Ki T 2(m-2) + K2 + 



PK 3 



j"2 j^o+i) 

Eliminating p(i) between Equations (21) and (27), we have 

1 



(27) 



A 



1 + 7-£o(™ + 2) 



(m+a-l)(l-m) T 2 ("+!) a(l-m)T 2 



Ki T 2(m-2) + ^2 + /? ^3 



J"2 y2(a+l) 



I + T 2 + T 2(a+1) 



(28) 

From Equations (26) and (28), we observe that when a > and m < 2, 
the positive cosmological constant is a decreasing function of time and 
approaches a small value in the present epoch. 
Some Physical Aspects of the Models: 

With regard to the kinematical properties of the velocity vector v l in 
the metric (19), a straight forward calculation leads to the following 
expressions for the scalar of expansion (9) and for the shear (a) of the 
fluid. 



(m + 2) 



\ 



y2(m-2) 



+ 



1 



(m + a - 1)(1 - m) T 2<Q+1 ) a(l - m)T 2 



(29) 
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T 2(m-2) p i 

(m + a- 1)(1 -m) + T 2 ( a+1 ) ~ a(l - m)T 2 

For a > and m < 2, the expansion factor is a decreasing function of 
T and approaches, asymptotically to zero with p and p also approaching 
to zero as T — > oo. 



1 



4. Particular Models 



If we set m = 2, then the geometry of the space-time (19) reduces to 
the form 



+ 



1 



dT 2 



T 2(3+64rf) 4(1 + 167^) ' (3 + 64tt^) 

+ T 4 dx 2 + T 2 dy 2 + (T 2 sin 2 y + T 4 cos 2 y)dz 2 -2T 4 cosy dx dz, (31) 
where /? is an integrating constant. 

The pressure and density of the model (31) are given by 

0(15 + 512tt£) 4(3 + 16tt£) (9 + 8ir£) 



8np - 

+ (320 
8vrp = 



+ 



3T 8(i+i6rf) 3(3 + 647r £) T 2 6 ( X + 167r£ ) 



^8(1+16tt£) 
5/3 



+ 



+ 



(3 + Utt£) T 2 4(1 + 16tt£) 
(1 + 8ir£) (3 + 8ir£) 



A 



T 8(i+i6rf) ( 3 + 64 ^) T 2 2(1 + 16tt^) 



+ A. 



(32) 
(33) 



4.1. Model I: 



Solution for £ = £ 



When n = 0, Equation (23) reduces to £ = £o = constant. Hence in 
this case Equation (32), with the use of (33) and (22), leads to 



8vr(l + 7 )p 



+ (32tt£ 



2/3(15 + 256n£) 4(15 + 32^£) 



+ 



Wit£ 



3T 8(i+i6rf) 3(3 + 647r£)T 2 3(1 + 16^) 



+ 



(34) 



_ T 8(i+i6rf) (3 + 64^)T 2 4(1 + 16tt£). 
Eliminating p(t) between (33) and (34), we obtain 

/3(15 + 512vrf - 15 7 ) 4 (3 + 16tt£ + 6(3 + 8^)7) (9 - 8tt£ + 3(3 + 8^)7) 



(1+ 7 )A 



3J , 8(1+16tt£) 



3(3 + 64n£)T 2 



6(1 + 1Gtt£) 
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+ (32tt£ ) 



4.2. Model II: 



P 



+ 



T 8(l+167rf) (3 + 647r £) T 2 4 (! + 167r ^) 



Solution for £ = 



(35) 



When n = 1, Equation (23) reduces to £ = £oP- Hence in this case 
Equation (32), with the use of (33) and (22), leads to 

1 



8vr 



1 + 7-4^0 



JL 



T 8(l+16^) + (3+64tt^) T' 2 4(1+16?^) 

2/3(15 + 256^£) 4(15 + 32tt£) Wtt£ 
3T 8(i+i6rf) + 3(3 + 64tt£)T 2 ~ 3(1 + 16tt£) 
Eliminating p(t) between (33) and (36), we obtain 

A = 



(36) 



8 



T 8(i+16tt<?) + (3+64vr£) T' 1 

4(15 + 32vr£) 



+ 



4(1+16tt^) 

16tt£ 



1 + 7-4^0 

2/3(15 + 256ttI) 

3T 8(i+i6rf) 1 3(3 + 64vr^)T 2 3(1 + 16^) 

5/3 8(l + 8vr£) (3 + 8tt£) 



T 8(i+i6rf) ( 3 + 64 ^) T 2 1 2(1 + 16t^) ' (37) 

From Equations (35) and (37), we observe that the positive cosmolog- 
ical constant is a decreasing function of time and approaches a small 
value in the present epoch. 
Some Physical Aspects of the Models: 

The expansion(#) and the shear(cr) in the model (31) are given by 



P 



+ 



1 



1 



T 8(i+i67rf) ( 3 + 64vr £) T 2 4 (! + 16ir £) 



a 



P 



+ 



(38) 



(39) 



_ T 8(i+i6rf) 1 (3 + 64^)T 2 4(1 + 16tt£). 

The expansion factor 6 in the model is a decreasing function of T. Since 
HmT^oof 0) hence the models do not approach isotropy for large 
values of T. The model have singularity at T = which is real physical 
singularity. 
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5. Special Models 

If we set m = 2 and t = Equation (18) leads to 

y/2 B (IB 

■ , = dt, 40 

v/(25 2 - £ 4 + 2/3) ' V ' 

which on integration gives 

B 2 = 1 + M sin(v / 2 t + 2N), (41) 

where M = a/2/3 + 1 and iV ia a constant of integration. Hence, we 
obtain 

A = B 2 = 1 + M sin(v / 2 t + 2N), (42) 
Using the transformations 

V2t + 2N = T, 
x = X, 

y = Y, 

z = Z, (43) 

the metric (1) takes the form 
dT 2 

ds 2 = — + (1 + M sinT) 2 dX 2 + (1 + M smT)dY 2 + 

[(1 + M sinT) sin 2 Y + (1 + M sinT) 2 cos 2 Y] dZ 2 

-2(1 + M sinT) 2 cosF dX dZ. (44) 
The pressure and density for the model (44) are given by 

_ 16\/2 vr£M cos T [35M 2 sin 2 T + 30M sin T - 2M 2 - 3] 

7TP ~ (1 + MsinT) + 12(1 + M sinT) 2 ' 

(45) 

[3 + 10M 2 + 2MsinT- HM 2 sin 2 Tl , . 

8vrp = 77 — ^-5 + A. 46 

F 4(1 + M sinT) 2 v ' 

5.1. Model I: Solution for £ = £ 

When n = 0, Equation (23) reduces to £ = £o = constant. Hence in 
this case Equation (45), with the use of (46) and (22), leads to 

_16\/2 7r( Mco8T [M 2 sin 2 T + 18M sinT + 14M 2 + 3] 
»n +7)P- ( 1 + MsinT ) + 6(l + MsinT) 2 ' 

(47) 
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Eliminating p(t) between (46) and (47), we obtain 



(1 + 7)A- l ^^ McosT [35M 2 sin 2 T + 30MsinT-2M 2 -3] 



(1 + MsinT) ' 12(1 + M sin T) 2 

[3 + 10M 2 + 2MsinT- 11M 2 sin 2 T] 7 
4(1 + MsinT) 2 ' 



(48) 



5.2. Model II: 



Solution for £ = £ p 



When n = 1, Equation (23) reduces to £ = £oP- Hence in this case 
Equation (45), with the use of (46) and (22), leads to 



8vrp 



1 , _ 2V2 goMcosT 
1 "•" ' (1+MsinT) 



[M 2 sin 2 T + 18M sin T + 14M 2 + 3] 
6(l + MsinT) 2 ' 

Eliminating p(t) between (46) and (49), we obtain 

1 



(49) 



8vrp 



1 + 7' 



cos 1 
(1+MsinT) 



[M 2 sin 2 T + 18M sin T + 14M 2 + 3] 
6(1 + MsinT) 2 

[3 + 10M 2 + 2MsinT- HM 2 sin 2 T] 

4(l + MsinT) 2 ' (50) 

From Equations (48) and (50), we observe that the positive cosmolog- 
ical constant is a decreasing function of time and approaches a small 
value in the present epoch. 
Some Physical aspects of the Models: 

The expansion(6*) and the shear(cr) in the model (44) are given by 



2^/2 McosT 
(1 + MsinT) 

I McosT 



(51) 



(52) 



^3 (1 + MsinT) 

When T -> then 6> -> 2\/2 M and when T -> f then 6> -> 0. Thus the 
expansion in the model starts at T = and it stops at T = ^. Since 
lim,T->oo7) 7^ 0, hence the models do not approach isotropy for large 
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values of T. The model has singularity at T = which is real physical 
singularity. 

6. Conclusions 

We have obtained a new class of Bianchi type IX anisotropic cosmolog- 
ical models with a viscous fluid as the source of matter. Generally, the 
models are expanding, shearing and non-rotating. In all these models, 
we observe that they do not approach isotropy for large values of time 
T. 

The cosmological constant in all models given in sections 3.1 and 3.2 are 
decreasing function of time and they all approach a small positive value 
as time increases (i.e., the present epoch). The values of cosmological 
"constant" for these models are found to be small and positive which are 
supported by the results from recent supernova observations recently 
obtained by the High - z Supernova Team and Supernova Cosmological 
Project ( Garnavich et al., 1998 ; Perlmutter et al, 1997, 1998, 1999; 
Riess et al, 1998; Schmidt et al. , 1998. Thus, with our approach, we 
obtain a physically relevant decay law for the cosmological constant 
unlike other investigators where adhoc assumption for the variation 
were used to arrive at a mathematical expressions for the decaying 
vacuum energy. 
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